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Abstract 

We obtain two new algorithms for partial fraction decompositions; the first is over al- 
gebraically closed fields, and the second is over general fields. These algorithms takes 
O(M^) time, where M is the degree of the denominator of the rational function. The 
new algorithms use less storage space, and are suitable for parallel programming. 
We also discuss full partial fraction decompositions. 
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1 Introduction 



The partial fraction decomposition of a one-variable rational function is very 
useful in mathematics. For example, it is crucial to obtain the partial fraction 
decon ipositiori of a r ational function in order to integrate it. Kovacic's algo- 
rithm lKovacic ( IQSfit ) for solving the differential equation y"{x) + r{x)y{x) = 0, 



where r{x) is rational, requires the full partial fraction expansion of r{x) over 
the complex numbers. 



The classical algorithm for partial fraction expansion relies on the following 
theorem. To make it simple, we consider rational functions in C{t). 

Theorem 1.1 // ai, . . . , a„ are n distinct complex numbers, mi, . . . , m„ are 

positive integers, and the degree of p{t) is less than mi + ■ ■ ■ + m„, then there 



Email address: maxima@brandeis.edu (Guoce Xin). 



Preprint submitted to Elsevier Science 



1 February 2008 



are unique complex numbers Aij, where 1 < i < n and 1 < j < rrii, such that 

'^'^ ttjA^Vr (1-1) 



The classical algorithm multiplies both sides by the denominator, and then 
equates coefficients to solve a large system of hnear equations for the v4j .,'s. 

The key observation for our first algorithm is that linear transformations will 
preserve the structure of a partial fraction expansion. This reduces the problem 
of finding all Aij to finding only Aij and assuming that Oi = 0. For finding 
Aij, we use the unique Laurent series expansion at t = 0. 

Denote by F{t) the left side of equation (1.1). Let M be the degree of the 
denominator of F{t), which is mi + m2 + ■ ■ ■ + mk. Compared with the classical 
algorithm for obtaining the partial fraction decomposition of F{t), our new 
algorithm has three advantages. This comparison is under the assumption of 
fast multiplication of (usually rational) numbers. In the following, when we 
say that an algorithm takes 0{M) time, we mean that the algorithm will do 
0{M) multiplications. 

(1) The new algorithm is fast. The classical algorithm needs to solve M linear 
equations in M unknowns, which take s VL{M ^) time using the Gaussian 
elimination algorithm. See (jSedgewick . 19881 p. 540, Property 37.1). But 



our algorithm only takes O(M^) time. 

(2) The new algorithm requires little storage space. The classical algorithm 
needs to record all of the coefficients in these M linear equations. 
So the storage space is about f2(M^). But our new algorithm needs only 
to record two polynomials of degree m, where m is the maximum of the 
mj's. So the storage space is at most 0(M). 

(3) The new algorithm computes the partial fraction expansion at different 
Oj's separately, so it is more suitable for parallel programming. 

For partial fraction decompositions in a general field we also have a fast 
algorithm by working in some quotient rings. The new algorithm applies to 
finding the full partial fraction decompositions and to evaluating generalized 
Dedekind sums. The theory gives an efficient algorithm for MacMahon's par- 



tition analysis (jXinl . 12004 Ch. 2.5). 



2 Partial Fraction Decompositions Over a General Field 



Let K be any field, and let t be a variable. It is well-known that the ring of 
polynomials K\t\ has many nice properties. Here we use the fact that K\t\ is 
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a unique factorization domain. We will use a quotient ring to derive a formula 
for partial fraction decomposition, which is the basis of our new algorithms. 



In what follows, the degree of an element r e K\t\, denoted by deg(r), is the 
degree of r as a polynomial in t. The degree of the polynomial is treated as 
— oo. We start with the division theorem in K\t\. 

Theorem 2.2 (Division Theorem) Let D,N e K[t\ and suppose D 0. 

There is a unique pair (p, r) such that p,r & K[t], N — Dp + r, and degr < 
degD. 

The r above is called the remainder of when divided by D. 

A rational function N/D with N,D E K[t] is said to be proper if degA^ < 
degD. A proper rational function is simply called a proper fraction. The unit 
1 is not proper, but is considered to be proper. It is clear that the sum of 
proper fractions is a proper fraction, and the product of proper fractions is a 
proper fraction. But the set of all proper fractions does not form a ring, for 1 
does not belong to it. 

By the division theorem, any rational function N/ D can be uniquely written as 
the sum of a polynomial and a proper fraction. Such a decomposition is called 
a ppfraction (short for polynomial and proper fraction) of N/D. UN — Dp+r 
with deg(r) < deg{D), then N/D = p + r/D is a. ppfraction. We denote by 
Poly(A^/L>) the polynomial part oiN/D, and by frac{N/D) the fractional part 
of N/D. 

Recall the following well-known result in algebra. 

Lemma 2.3 Let N, D E K[t] with D 0. If D = Di ■ ■ ■ D^ is a factorization 
of D in K[t], and all the Di are pairwise relatively prime, then N/D can he 
uniquely written as 

N ri rk , ^ 

where ri is a polynomial of degree smaller than deg{Di) for all i, and p equals 
the polynomial part of N/D. 

Wc call the above decomposition the ppfraction expansion of N/ D with respect 
to {Dij . . . , Dk). To find such decomposition, it suffices to find p and ri, . . . , r^. 
It is easy to find p. In finding r^, it is convenient to use the concept of quotient 
ring K[t]/{Di), where (A) is the ideal generated by Di. 

Recall that D' + (D) has a multiplicative inverse in K[t]/ (D) if and only if D' 
is relatively prime to D. Moreover, if ~ 13 ~^ jyy which is not necessarily 
a ppfraction expansion, then 1/D' + (D) = r + {D). The last fact we will use 
is that for any polynomials A^i and N2, Ni + (D) — N2 + (D) if and only if 
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Frac(A^i/D) = Frac(A^2/i^). 

Suppose that D = DiD' and that Di and D' are relatively prime. Then we 
have a ppfraction of N/D with respect to {Di, D'): 

N/D = Po\y{N/D) + n/Di + r'/D'. 

In such a decomposition, we call ri/Di the fractional part of N/D with respect 
to Di, and denote it by frac{N/D, Di). If -Di = {t — a)"* for some a E K, then 
we simply denote it by frac{N/ D,t — a). Clearly Frac(A^/D, 1) is always 0, 
and Frac(A^/D, Di) is always a proper fraction with denominator Di. Also we 
have the following relation: 

n + {D,) = ^ + {D,). 

Thus to find ri, we pick a representative of N/D' + (-Di), and then find its 
remainder when divided by Di. 

Theorem 2.4 For any N,D E K[t] with D ^ 0, if Dx,...D^ e K\t\ are 
pairwise relatively prime, and D = Di - ■ ■ Df;, then 



N ^ , fN\ ^ fN ^\ ^ fN ^ 

- = Poly i^—j + Frac (^-, D^j+--- + Frac (^-, D^ 



is the ppfraction expansion of N/D with respect to {Di, . . . , Dk). Moreover, if 
l/(L>iA) = Si/Di +Pi/Di, then 

Frac{N/D, D^) = Frac(iVs2S3 • • • Sfc/^i)- 



PROOF. For the first part, suppose that 

N ri ru , , 

is the ppfraction expansion of N/D with respect to (Di, . . . , Dy/). Let D' — 
D2 • • • Dk- Then Di and D' are relatively prime and r-^/D^-^- • ■+rk/Dk = r'/D' 
is a proper fraction with denominator D'. By the uniqueness of ppfraction of 
N/D with respect to {Di,D'), we have ri/Di = frac{N/D,Di). Similarly 
ri/Di = Frac(iV/L), A) for all i. 

For the second part, multiplying both sides of equation (2.2) by D, and think- 
ing of this as an identity in the quotient ring K[t]/{Di), we get 

TV + (D,) = pD,D' + nD^ ■■■Dk + --- + rkD, ■ ■ ■ D^-i + {D^} = r^D' + {D^). 
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Now multiply both sides of the above equation by 1/D' + (Di), we get 

N 

ri + (Di) = — + {Di) = NS2S3 ■■■Sk+ (D,). 

U2 - ■ ■ JJk 

Therefore 

Frac -— = Frac 



Theorem 2.4 is the basis of our new algorithms. Let M be the degree of the 
denominator of a rational function. We will give an O(M^) algorithm for 
finding the partial fraction decomposition based on the above theorem. 

li D = ap^^ ■ ■ ■p^'', where a E K, is a factorization of D into monic primes 
in K[t]j then . . . jP™* are pairwise relatively prime. Let Di = p^' , and let 
Tj be a polynomial with deg(rj) < deg(Dj). Then every Vi/Di can be uniquely 
written in the form J2Y=i ^j/Pi with deg{Aj) < deg(pi) for all j. The partial 
fraction expansion of N/ D is the result of applying the above decomposition 
to the ppfraction of N/D with respect to {Di, . . . , Dk). In this case, we can 
use the following lemma to reduce the problem to computing only the partial 
fraction expansion of l/{piPj) for all i ^ j. 

Lemma 2.5 Let p,q E K[t] be relatively prime polynomials. If r and s are 
two polynomials such that 1/ipq) — r/p + s/q, then for any positive integers 
m, n, 

= - E f"^ ^ '] r-^^ + - E ^ r^s-q^. (2.3) 



PROOF. Using the formula 1/ipq) = r/p + s/q, we have 

111 r s 

+ 



pTnqTi pq pTn—lqU—l pmqn—l pTn—lqu' 

If we let A{m,n) — l/(p"*g"), then the above equation is equivalent to 

A{m, n) = rA{m, n — 1) + sA{m — l,n). 

Using this recursive relation, we can express A{m,n) in terms of A{0,j) and 
A{i, 0), where 1 < j < n and 1 < i < m. 

Either using induction or a combinatorial argument, we can easily get 
A{m, n) VA(m - i, 0) + E ("^ ^ ^Vh^^AiO, n-j). 

i=Q \ ^ / j=0 \ J J 

Equation (2.3) is just a restatment of the above equation. 
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3 Pcirtial Fraction Decompositions in C{t) 

In this section, K is an algebraically closed field (e.g., the field of complex 
numbers C). Partial fraction decomposition in this situation is simple, since 
every polynomial in K[t] can be written as a product of linear factors t — a 
for a & K. 

The key idea to our new algorithm is that linear transformation will not change 
the structure of a partial fraction decomposition. This can be illustrated by 
the following example. 

The partial fraction expansion of /(t) is A/{t — a) + B/{t — h) if and only if 
the partial fraction expansion of /(i + c) is A/ + c — a) + S/ + c — 6) . So we 
can compute the partial fraction expansion of f{t + a), and after that, replace 
t with t — a. 

Let h & K and let rf, be the transformation defined by f{t) = f{t + h) for any 
fit) e K\t\ or fit) e Kit). Then T{, is clearly an automorphism of K\t\ and 
of K{t), and its inverse is r_6. The following properties can be easily checked 
for any p, g G K[t] and b E K. 

(1) p is prime in K[t] if and only if p is. 

(2) n gcd{p,q) = gcdin p,n q)- 

(3) deg(rb p) = dcg(p). 

(4) p/q is a proper fraction if and only if T{, p/q is. 

Thus for any N,D E K[t] with D ^0, N/D = p + n/Di H h Vk/Dk is the 

ppfr action expansion of N/ D if and only if N/ D = {u p) + {Tb ri/ Di) + ■ ■ ■ + 
i^b fkl Dk) is a ppfraction expansion of r?, N/D. The partial fraction expansion 
can be obtained by first computing the partial fraction expansion of Tb N/D, 
then applying r^b to the result. Choosing b appropriately can simplify the 
computation. The above argument gives us the following lemma. 

Lemma 3.6 For any N, D, Di e K[t] with D ^Q, D/D^E K[t\, and 
gcd{Di, D / Di) — 1, we have 

frac{N/D,Di) = T-b fracin N/D^n Di). 

Let [t™] be the map from to K\t\ given by replacing with for all 

n >m. More precisely, 

m—l 

n>0 n=0 
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where ai G K for all i. The following properties can be easily checked for all 
f,9^K[[t]]- 

(1) {f + g)= ft™] / + m g. 

(2) \mfg)= \m\nf\ng). 

(3) If < A; < m then ft™] t'^f = [t™"^] /. 

(4) If ^7(0) ^ 0, then [t^] f/g = [t-] ( [t™] // [t™] 

The main formula for our algorithm is the following, which is a consequence 
of Theorem 2.4. But we would like to prove this result by using Laurent series 
expansion. 

Theorem 3.7 Let N,D e K[t] and let D = f^E with E e K[t] and E{Q) ^ 0. 
Then 

N{t) 



PROOF. Since E{0) ^ 0, and E are relatively prime. Let 

^-v(t) + '-^ + ^ (31) 

be the ppfraction of N/D with respect to (t"^,E). Thus deg(r(t)) < m, and 
r{t) = t"'frac{N/D,t"'). 

Because K{t) can be embedded into the field of Laurent series K{{t)), equation 
(3.1) is also true as an identity in K{(t)). On the right-hand side of equation 
(3.1), when expanded as Laurent series in K{{t)), the the second term contains 
only negative powers in t, and the other terms contain only nonnegative powers 
in t. Therefore, r{t)/t"^ equals the negative part of N/D when expanded as a 
Laurent series. More precisely, for i = 1, . . . , m, we have 

This is equivalent to [t"^-']N{t)/ E{t) = [t™~^]r(t) for i = 1, . . . , m. Now r{t) 
is a polynomial of degree at most m — 1, and N(t)/E(t) G -ft'iit]], so 

r(t) ^ in If. 

Remark 3.8 The idea of using Laurent series expa nsion to obta in part of 
the partial fraction expansion appeared in the proof of i Gessel 1991 . Theorem 

44). 
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Gessel observed that this same idea can also be used to compute the polyno- 
mial part of a rational function, and that it is fast when the polynomial part 
has small degree. 

Proposition 3.9 If R{t) is a rational function in K{t), then the polynomial 
part P{t) can he computed by the following equation. 

t-ip(t-i) = frac{t-^ R{t-^),t = 0). 



PROOF. Let R{t) = P{t) + N{t)/D{t) be the ppfraction of R{t), and let 
p — deg(P), d — deg{D), and n — deg{N). Then n < d. Now we have 

t-^R{t-') = t-^P{t-^)+t-^N{t-^)/D{t-^) = r^P(r^) + t'^-"-^7V(t)/(L>(t)), 

where D{t) = t'^D{t-^), and similarly for N{t). 

Apply ppfraction expansion to the second term. Since D{t) has nonzero con- 
stant term, it is relatively prime to f^^^. Now it is clear that t~^P{t~^) is the 
fractional part of t~^R{t~^) with respect to t'^'^^. 

Example 3.10 It is easy to check that 

, t3 + 2t2-3t + 4 ^ -8 + 19t 



^2-4^ + 2 t'^-At + 2' 

Now we compute the polynomial part of R{t) by Proposition 3.9. 



l + 2t-3i^+4i3 
^- t^l-At + 2t^) 
t'Frac{t-'R{t-'),t^)= p^2^ 1 + 2^-3*^ + 4^3 



(1-4^ + 2*2) 



21 1 + 2t 



So we obtain that the polynomial part of R{t) is i -|- 6. 



Note that when expanded as Laurent series in we have 

1 ™o~^ /m - — ] + 'i\ 



(t-a,)". ^„ \ j laT*' 

Hence by Theorem 3.7, we get 
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Corollary 3.11 Let N G K[t] and D = t'^'^it - ai)""^ ■ ■ ■ {t - 0^)"^'= with all 
the ai's distinct and not equal to 0. Then 

r«Frac (—, = Nsi--- s^, 



where 



j=0 



■mi+j ' 



Therefore, combining Theorem 2.4, Lemma 3.6 and Corollary 3.11, we obtain 
an algorithm for computing the partial fraction decomposition of a proper 
rational function of the general form 



(t-ai)"i---(t-afe)-'=' 

(1) Let 5* := 

(2) For i from 1 to A; do G{t) := F{t + ai), S := S + r_„Trac(G'(t), t'^O next 
i. 

(3) Return S. 

It was stated in (|xinl . 12004 Ch. 2.4) that the computation of Frac(G(t), t™') for 
all i will take time 0{k{m\'^^+m2^^ + - ■ ■+m^'^'^)). However, this estimate does 
not show O(M^) time. For example, M = 2k, and mi = 1712 = ■ ■ ■ = rrik-i, 
and rrik = M — k + 1. 

Theorem 3.12 Let M be the degree of the denominator of a rational func- 
tion. The above algorithm for partial fraction decomposition can be executed 
in O(M^) time. 

The proof of this theorem, which will be given later, uses the fact that ma- 
nipulations in K[t]/ (t™) are fast. Now let us estimate the computational time 
of manipulations in the quotient ring K[t]/{t'^). 

The follo wing is a well-kno wn result by the method of divide and conquer. 
See, e.g., (Sedgewick, 19881 Property 36.1). 

Proposition 3.13 Let R{m) be the time for computing the product of two 
polynomials of degree less than m. Then R{m) = 0{m^'^^). 

Remark 3.14 In the proof of Theorem 3.12, we only need the obvious upper 
bound R{rn) = Oim?). The above proposition shows that our algorithm can be 
accelerated. 

Since most of our estimations use the method of divide and conquer, it is 
better to introduce it here. In what follows, we shall always assume that m is 
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a power of 2 for simplicity. The estimation of R{m) follows from the following 
observation. Bisect P{t) as P{t) = Pi{t)+t'^/'^P2{t), and bisect Q{t) as Q{t) = 
Qiit)+f"'/^Q2it). Then 

PQ = p,Q, + r/'((Pi + Q,)(P2 + Q2) - PiQi - P2Q2) + rP2g2. (3.2) 

which shows that we need only three polynomial multiplications. This gives 
the recurrence R{m) = SR{m/2) and that R{m) = 0(m'°S2^) = 0{m^-^^). 

Lemma 3.15 Let P(t) and Q{t) he two polynomials of degree m — 1. Then 
ft™] P{t)Q{t) can he computed in no more than R{m) time. 

The proof of this lemma is trivial. 

Lemma 3.16 The computation of [t™] P/Q, where Q{Q) ^ 0, takes no more 
than 2R{m) time. 



PROOF. We use the method of divide and conquer. Let T(m) be the time 
for the computation in question. 

Suppose that [t'"] P/Q = Z. Bisect P,Q,Z as P ^ Pi + t"'/^P2, Q ^ Qi + 
t"'/^Q2, and Z = Zi + t'^/^Z2. Then 



Obviously [W^^ Pi/Qi = Zi. So it will take T(m/2) time to find Z^. To find 
Z2, we use the formula (from direct algebraic computation). 

Therefore, we get the recurrence T(m) = T(m/2) +2i?(m/2) + T(m/2), where 
the first summand is for Zi, and the rest is for Z2. Using this recurrence, it is 
easy to see that T(m) is no more than 2R{m). 



To apply the above lemma, we need the expanded representation of P{t) and 
Q{t). 

Lemma 3.17 Suppose that Q(t) is the product of m linear factors. Then it 
takes no more than R{m) time to expand Q{t). 



PROOF. Let U (m) denote the time for expanding the products of m linear 
factors. Factor Q{t) as Q{t) = Qi{t)Q2{t), where Qi{t) consists of the first 
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m/2 factors. Then it will take U{m/2) time to expand Qi{t), and U (171/2) 
time to expand Q2{t), and then R{m/2) time to get the final expansion. Thus 
U{m) = 2U{m/2) + R{m/2). This implies that U{m) is approximately equal 
to R{m). 

We might be able to speed up the expansion in the above lemma by the 
following lemma. 

Lemma 3.18 The expansion of (t — a)"* takes 0(m) time. 
PROOF. This lemma follows from the binomial theorem 



and the fact that the ratios of consecutive summands are simple. 

Using the binomial theorem, it is easy to see the following. 

Lemma 3.19 Suppose the degree of N{t) is less than M. Then the expansion 
of [t™] N(t+a) takes no more than 0{m-the number of nonzero terms in N{t)), 
which is no more than 0{Mm) time. 

Now we estimate the computational time for frac{N/D, mo) in Theorem 3.7. 
Proposition 3.20 Let M = niQ + mi + ■ ■ ■ + m^. Then it takes 0{Mm^-^^) 



PROOF. We first expand the denominator by grouping every (mg — 1) factors 
together. So we have about (M — mo)/(mo — 1) groups. It takes R{mo) time 
for expanding the products for each group, and then about (M — mo) / (mo — 1) 
multiphcations when taking \t'^°'\ . Thus the total time for this expansion is 
about 2M/moR{mo). 

Denote by E(t) the resulting expansion. Now it will take about 2i?(mo) time 
to compute [t™"] Therefore the total time for the final answer is 

0{M/mo)R{mo) = 0(Mm[]-5«). 

Proof of Theorem 3.12. It will take 0{Mmi) time for finding the expansion 
oiN{t + ai), and wiU take 0{Mmf^^) time for finding ft™^] t''''^N{t + ai)/D{t + 
Oj). Therefore, the total time for finding frac{N/D,x — ai) takes 0{Mmi) 




time to compute [t™""] 



mi 



{t-ai)"'i-it-ak)"'k ■ 
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time. Summing on all i, we see that it takes 0{M^) time to find the partial 
fraction decomposition of N/D. 



This new algorithm also enables us to work with some difficult rational func- 
tions by hand. 

Example 3.21 Compute the partial fraction expansion of fit), where 

t 



fit) 



{t+l)\t-l)^{t-2f 



Solution. Clearly, the polynomial part of fit) is 0. Although applying Corol- 
lary 3.11 is faster, we compute the fractional part of /(t) at t = — 1 and 
t = 1 differently. For the fractional part of f{t) at t = —1, we apply r_i, and 
compute Frac(/(t — l),t^) by Theorem 3.7. We have 

evrzc{f{t-i),e)= 



(t-2)3(t-3)5 



-8 + 12t)((-3)5 + 34-5t) 



8-35(l-19/6t) 

^ 2n (t-l)(l + 19/6t) ^ I m 

' ' 8-3^ 8-35^ Q 

Thus 

Frac(/W, {t + m = - ,3.34 + 1)2 " ^TTWTTy 
Similarly, we can compute the fractional part oi f{t) a,t t — 1. We have 
t'frac{f{t + l),t')^ \t'] ^ + ^ 



(t + 2)2(t-l)5 



31 ^ + 1 



Thus 

Frac(/(i),(i-lf) 



{t^ + t + A){-10t^ + 5t - 1) 

H^l ill 

' ' -4 + 16t - 21t2 

-^\t'] {t+i)ii + At-^t' + m') 

__(l + 5t + -t^). 



1 5 59 



4(t-l)3 4(t-l)2 16(t-l)' 
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The fractional part of f{t) at t = 2 can be obtained similarly, but it is better 
to use Corollary 3.11. In fact, this computation becomes quite complicated. 
Although it is still possible to work by hand, we did use Maple. 



t^Frac{f{t + 2),t^) 



729 



(l - 3t + 6t^ - + 15t^) 




2 19 13 2 593 o 2689 



Applying Theorem 2.4, we get the partial fraction expansion of f{t), which is 
too lengthy to be worth giving here. 



4 Algorithm for a General Field and Full Pcirtial Fraction Decom- 
positions 

When 7^ is a general field, e.g., the field of rational numbers Q, linear trans- 
formations will not help. Manipulations in K[t]/ {Di{t)) are not as good as the 
case of Di{t) = x"\ But we still have an O(M^) algorithm. 

Proposition 4.22 Suppose that Di,...,Dk G K[t] are pairwise relatively 
prime, and D = Di- ■ ■ D^. If deg{N) < deg{D), then the ppfraction decom- 
position of N/D with respect to D^, . . . ,Dk can be computed in 0(deg(D)^) 
time. 

The proof of this proposition will be given later. Now assume that = p"% 
and deg(Dj) = mj. It is easy to show that the partial fraction decomposition 
of Ti/Di, where deg(ri) < rrii, can be computed in 0{m1) = 0{Mmi) time. 
Thus the above argument and Proposition 4.22 will give us the following. 

Theorem 4.23 Suppose that deg(iV) < deg{D), and we are given a factor- 
ization D = pi^ ■ ■ ■pi'' of D into primes in K[t]. Then the partial fraction 
decomposition of N/D takes 0(deg(D)^) time. 

In order to prove Proposition 4.22, we need to estimate manipulations in 
K\t\/ {D) for a given polynomial D. In most situations, we need the unique 
representative of + {D) that has degree less than degD. We denote by 
\D~\ N this representative, which is also known as the remainder of A^ when 
divided by D. 

The following estimations are obvious. The computational time refers to the 
number of multiplications of two elements in K. Time spent on additions is 
omitted. 
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(1) The computational time for expanding PQ for any two polynomials P 
and Q is no more than (deg(P) + l){deg{Q) + 1) time. 

(2) For two given polynomials and D, the division algorithm for finding p 
and r such that = pD + r with deg(r) < deg(i5) takes no more than 
(deg(Ar) - deg(L>)) ■ (deg(L>) + 1) time. 

(3) Suppose that P and Q are two polynomials of degree less than deg(L'). 
Then the computation of [Z^] PQ takes no more than 2deg(D)^ time. 

Lemma 4.24 Suppose P is relatively prime to D and deg(P) < deg(Z^). 
Then the computation of \D'\ 1/P takes 0(deg(D)^) time. 

This estimat ion is ob tained by the extended Euclidean algorithm for polynomi- 



als. See, e.g.. lMoenck (1973,), which says that an 0(deg(D) \o^{D)) algorithm 



exists, where r is a fixed number. 

Lemma 4.25 Suppose that deg(Di) = m and deg(i52 ■ ■ ■ -Dfc) = M . Then the 
computation of \Di \ D2 - ■ ■ Dk takes no more than 4(M + m)m time. 



PROOF. Denote by V{M) the computational time described in the lemma. 
We shall prove that V{M) < max{ 4Mm — 2m^, Mm + }, which implies 
the lemma. The proof is in two parts. The first part deals with the case when 
M < 2m, and the second part deals with the case when M > m. Note that 
there is an overlap. 

We first show that the expansion of Z^2 ■ ■ ■ -Dfc takes no more than M^/2 time 
by induction on M. This claim is clearly true for small M, e.g., M = 1,2. Now 
suppose the claim is true for all / < M. Then the expansion of - ■ ■ can 
be obtained by first expanding D2 ■ ■ ■ -D/t-i (of degree Mi), then multiplying 
it by Dk (of degree M2). The computational time is (by induction) no more 
than Ml/2 + M1M2 < (Mi + M2)V2 = MV2. Thus for M < 2m, we get 
V{M) < M^/2 + m^ < Mm + m^ by expanding - ■ ■ D^, and then taking 

We claim that for all M > m, V{M) < AMm - 2m^, and prove the claim 
by induction on M. The claim follows from the inequality M^/2 + m^ < 
4Mm — 2m^ when M < 7m by the first part. For M > 7m, we can separate 
D2 - ■ ■ Dk into two products of degree Mi and M2 respectively. We can assume 
that Ml > m and M2 > m, for otherwise, the degree of one of Di is larger 
than 5m, in which case the claim is easily seen to be true. Now we compute 
the remainder of each product, and then compute the resulting product and 
compute the remainder. This process takes time 

V{Mi) + V{M2) + 2m^ < 4Mim - 2m2 + 4M2m - 2m^ + 2m^ = 4Mm - 2ml 
This completes the proof. 
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Proposition 4.26 Suppose that deg(D) = M , deg(A^) < M , deg(Di) = m, 
and that D = Di ■ ■ ■ is a factorization of D into relatively prime factors. 
Then the computation of \Di \ N/{D2 ■ ■ ■ Dk) takes 0{Mm) time. 



PROOF. Wc first compute \Di] D2 - ■ ■ Dk, and denote the result by D'. This 
step takes no more than 4Mm time. Then we compute \Di] N, and denote 
the result by N'. This step takes no more than Mm time. Finally we compute 
\Di] N'/D'. This step takes 0{m^) = 0{Mm) time by Lemma 4.24. So the 
total time is AMm + Mm + 0{Mm) = 0{Mm). 



Proof of Proposition 4.22. By Theorem 2.4, the numerator of Frac{N/ D, D^) 
is given by 

r, = \D,^ND,/{D^D2■■■Du). 
The computation of r, takes 0{M deg(Dj)) time by Proposition 4.26. Summing 
on all i we get the total computational time for the ppfraction of N/ D with 
respect to Di, . . . , D^, which is Y!1=i 0{Mmi) = 0{M^). 

Example 4.27 Compute the fractional part of f{t) with respect to —t + 2, 
where 

■^^^^ " (^2 -2t- 1)2(^2 -t + 2y 



Solution. Let p{t) = t'^-t+2. Then we need to compute \p{t)'\ t^/{f-t-lf. 
In the following computation, we shall always replace t"^ with t — 2. 



(i2-2i-l)2 '^^^'(-^-3)2 '^^^'7(t+l) 
\Pyt)\ _28 ' 

where we used the fact that {t + l){t — 2) — t'^ — t — 2 — p{t) — 4. Therefore 

— 9 

Frac(/(t),p(i)): 



-28(^2 - t + 2)' 



In Maple, the full partial fraction expansion of a rational function will involve 
a form like 

a=root of pit) J=l ^ ' 

where p{t) is a prime polynomial, and hj{t) will be a polynomial of degree 
less than deg(p(t)). This expansion is useful in some situations. We can get 
this kind of expansion by applying Theorem 3.7. This is best illustrated by an 
example. 
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Example 4.28 Compute the full partial fraction expansion of f{t), where 

t 



fit)- 



{t''-t-iy{t^-t + 2)' 



Solution. Suppose that a is a root of tlic prime polynomial p{t) := — t—1. 
Since K{a) is a field, and = a + 1, we can use this relation to get rid of all 
terms containing a"^ for n > 2. Because p{t) is a prime polynomial, a can only 
be a simple root of p{t). Then t divides p{t + a) and p{t + a)/t has nonzero 
constant term. In the present example, 

p(t + a) = {t + af - (t + a) -1 = t{t + 2a- 1). 

Note that the constant term of p{t + a) is always 0. 

Clearly, r^it"^ — t + 2) has nonzero constant term, for otherwise t'^ —t + 2 will 
not be relatively prime to p{t). In the present situation, 

{t + af - {t + a) + 2^t^ + {2a - l)t + 3. 

By Lemma 3.6 and Theorem 3.7, we can work in i^(Q;)[[t]]. 

\t'^t'f{t + a)= \t'] 



{t + 2a- l)2(t2 + (2a - l)t + 3) 

= — r^^i ^ + ^ 

is' ' 1 + (2a- l)llt/15 

= ^ \t^] {t + a)(l - ll(2a - l)t/15) 
lo 



1 -11a + 7 

= —a + ^ 

15 152 

Thus the fractional part of f{t) at a that satisfies p{a) — can be written as 

a _|_ ~ 11a) 



15(i-Q;)2 225(t-a)' 

Similarly, the fractional part of f{t) at P that satisfies P'^ — P + 2 — can be 
written as 

Together with the fact that the polynomial part of f{t) is clearly 0, the full 
partial fraction expansion of f{t) is hence 



m = E 



a _|_ ~ 11a) 



15(t - q;)2 225(t - a) 



-9^2=0 63(i - P) 
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Of course we can first compute the partial fraction decomposition of f(t) and 
then compute its full partial fraction decomposition. 
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